Magnetic Excitations in the Site-Centered Stripe Phase: 
Spin Wave Theory of Coupled Three-Leg Ladders 
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The success of models of coupled two-leg spin ladders in describing the magnetic excitation spec- 
trum of La2_a:Baa;Cu04 had been interpreted previously as evidence for bond-centered stripes. In 
a recent article, however, we have determined the magnetic coupling induced by the charge stripes 
between bond- or site-centered spin stripes modeled by two- or three-leg ladders, respectively. We 
found that only the site-centered models order. We further indicated excellent agreement of a 
fully consistent analysis of coupled three-leg ladders using a spin wave theory of bond with the 
experimental data. Here, we provide a full and detailed account of this analysis. 

PACS numbers: 74.72.-h, 74.20.Mn, 75.10.-b, 75.25.-j 
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I. INTRODUCTION 

Twenty years after the discovery, the mechanism of 
high-T(. superconductivity in the copper oxide materi- 
als is still considered one of the most important out- 
standing problem in contemporary physics^— . The ma- 
terials are described by mobile charge carriers (holes) 
doped into a quasi-twodimensional spin 1/2 antifer- 
romagnet^i^. Inelastic neutron scattering experiments 
have revealed a magnetic resonance peak^i^ and, in 
some compounds, periodic modulations in the spin and 
charge density (stripes Tranquada et ali^ found 
that the magnetic excitation spectrum of stripe or- 
dered Lai.875Bao.i25Cu04 looks similar to disordered 
YBa2Cu306-i-a~ or Bi2Sr2CaCu208+^, and observed 
that the data are consistent with bond-centered stripes 
modeled by two-leg ladders. This experiment is consid- 
ered of key importance for the field, as it it may provide 
the decisive hint as to within which framework copper 
oxide superconductors may be understood. 

With regard to such a framework, there is no consen- 
sus at present, but a fierce competition among differ- 
ent schools of thought. One of these school o^^d^d'^d^ 
attributes the unusual properties of the doped, two- 
dimensional antifcrromagncts to their propensity to form 
stripes, or their proximity to a quantum critical point 
(QCP) at which stripe order sets in. The resulting 
picture is highly appealing. Static stripes have been 
observed^ii^ only in certain compounds, most notably 
La2-2:Sr2:Cu04 at a hole doping concentration a; = i, 
and are known to suppress superconductivity. On the 
other hand, the mere existence of stripes would im- 
pose an effective one-dimensionality, and hence provide 
a framework to formulate fractionally quantized excita- 
tions. This one-dimensionality would be roughly consis- 
tent with an enormous body of experimental data on the 
cuprates, including the electron spectral functions seen 
in angle- resolved photo emission spectroscopy (ARPES). 
The charge carriers, the holons, would predominantly re- 
side in the charge stripes, as they could maximize their 
kinetic energy in these antiferromagnetically disordered 



regions. In the spin stripes, by contrast, the antiferro- 
magnetic exchange energy between the spins would be 
maximized, at the price of infringing on the mobility of 
the charge carriers. Most importantly, the spin stripes 
would impose a coupling between the charge stripes, 
which would yield an effective, pairwise confinement be- 
tween the low-energy spinon and holon excitations re- 
siding predominantly in the charge stripes. The mech- 
anism of confinement would be similar to that of cou- 
pled spin chains or spin laddersii"— . The holes would be 
described by spinon-holon bound states, and the domi- 
nant contribution to the magnetic response measured in 
Tranquada's as well as all other neutron scattering ex- 
periments would come from spinon-spinon bound states. 

The similarity of the "hour-glass" spectrum shown in 
Fig. 4b of Tranquada et ali^ (which is reproduced for 
comparison in Fig. [TJj) with the "elephants trousers" 
observed by Bourgcs et alj^ (which are reproduced for 
comparison in Fig. [Ud) provides the most striking ev- 
idence in favor of the picture advocated by this school, 
which attributes the anomalous properties of generic, dis- 
ordered CuO superconductors to the formation of dy- 
namic (rather than static) stripes, which fiuctuate on 
time scales which are slow compared to the energy scales 
of most experimental probes. This picture is considered 
to receive additional support by Xu et al^, who ob- 
served that the magnetic response of Lai.875Bao.i25Cu04 
at higher energies is independent of temperature, while 
the stripe order melts at about Tst ~ 54K. Measure- 
ments on 'untwinned' samples of YBa2Cu306.6j where 
one would expect the dynamical stripes to orient them- 
selves along one of the axis, however, exhibit a strong 
anisotropy in the response only at energies below the 
resonance^!, while the response is fourfold rotationally 
symmetric at higher energies^^. We believe, however, 
that this only indicates that the formation of stripe cor- 
relations, be it static or dynamic, is a low energy phe- 
nomenon, while the high energy response probes itiner- 
ant antifcrromagncts at length scales on which the stripes 
are essentially invisible. 

An extremely appealing feature of the experiment by 
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FIG. 1: (Color online) (a) Superpositions of cuts along {kx, n) 
and (tt, ky) for the lowest mode cj(fc) obtained with the bond 
operator spin wave theory of coupled three-leg ladders pre- 
sented here (red) superimposed with the experimental data 
obtained by inelastic neutron scattering by Tranquada et ali^ 
(black), (b) The neutron data as originally presented, with a 
triplon dispersion of a two-leg ladder superimposed (red line) 
(Reprinted by permission from Macmillian Publishers Ltd: 
Nature 429: 534-538, ©2005). 




FIG. 2: Overall momentum dependence of the magnetic 
response of superconducting YBa2Cu3 06.85 as reported by 
Bourges et al^. There is no stripe ordering in this com- 
pound (From Bourges et al, Science 288: 1234-1237 (2000). 
Reprinted by permission from AAAS). 



Tranquada et alt^ is that it immediately suggests a 
model of ferromagnetically coupled two-leg ladders, as 
the upper part of the measured spectrum agrees strik- 
ingly well with the triplon (or spinon-spinon bound state) 
mode of isolated two- leg ladders (see Fig.[T]D). The exper- 
iment hence appeared to point to bond-centered rather 
than site-centered stripes, (i.e., stripes as depicted in 
Fig. [3K rather than Fig. Wp)-, and thereby to resolve a 
long outstanding issue. 

This interpretation received support by theoretical 
studies^^"— . Vojta and Ulbricht^ used a bond operator 
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FIG. 3: Spin model for (a) bond-centered and (b) site- 
centered stripes in CuO superconductors. 



formalisni^^ similar to ours to study a spin-only model 
of stripes^I of coupled two-leg ladders (as depicted in 
Fig. [3^) with J|[ = J, took into account a bond-boson 
renormalization^^ of J, and assumed a value J' for the 
ferromagnetic coupling between the ladders which is large 
enough to close the spin gap of the ladders, i.e., to in- 
duce long range magnetic order. Within their approxi- 
mations, a value of J' = — 0.06J is sufficient. This value 
is not consistent with previous studies^^Ti^^ but as no 
method to calculate or even estimate the true J' induced 
by charge stripes had been available, it did not seem a 
problem at the time. The spectrum they obtained agrees 
well with experimental data measured by Tranquada et 
ali^, and hence appeared to justify their assumptions 
a posteriori. They concluded in favor of bond-centered 
stripes. This conclusion was independently strength- 
ened by Uhrig, Schmidt, and Griiningei^, who used the 
method of continuous unitary transformations to study a 
model of ferromagnetically coupled two-leg ladders, and 
observed that the critical value of J[. can be significantly 
reduced if a cyclic exchange term Jcyc on the ladders is in- 
cluded^. They likewise fine-tuned J' to the QCP where 
the gap closes and long-range magnetic order ensues, and 
reported good agreement with experiment. 

On the other hand, Seibold and Lorenzana^^^ calcu- 
lated the magnetic response for a range of dopings within 
the time-dependent Gutzwiller approximation, and found 
good agreement with the measured data for both bond- 
and site-centered stripe models. 

In a recent article^, we investigated whether it is rea- 
sonable to assume that the ferromagnetic coupling J' in- 
duced by the charge stripe between the spin stripes mod- 
elled by two-leg ladders is sufficiently to induce long range 
order. There are several estimates for the critical value 
required if the coupling between isotropic ladders is 
antifcrromagnetic in the literature. Gopalan, Rice, and 
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SigristSS find « 0.25J in a simple mean-field treat- 
ment of bond-bosons. Quantum Monte Carlo (QMC) 
calculations by Tworzydlo et alM yielded = 0.30(2) J, 
a value subsequently confirmed by Dalosto and Riera^. 
We redid the mean-field calculation of Gopalan et al^ 
for ferromagnetic (FM) couplings < 0, and found that 
within this approximation, the absolute value of is 
independent of the sign of the coupling. QMC calcu- 
lations by Dalosto et al^, however, indicate that the 
true value is at least Jc = — 0.4J (see Fig. 6b of their 
article). The physical reason why a significant coupling 
between the ladders is required to induce magnetic order 
is that the individual two-leg ladders possess a gap of or- 
der A « J/2. As a cyclic exchange term Jcyc ~ 0.25J 
reduces this gap by a factor of two2^, we expected that 
would likewise be reduced by a factor of two. We hence 
concluded that a FM coupling of at least somewhere be- 
tween = — 0.2J and — 0.4J is required, depending on 
the strength of a possible cyclic exchange term. 

The value we obtained for the ferromagnetic coupling 
induced by the charge stripes between the spin stripes 
through exact diagonalization of small clusters with and 
without charge stripes, however, is J' = —0.05 J—. The 
details of this calculation are given in Appendix |B] The 
coupling is hence insufficient to induce order in a model of 
coupled two-leg ladders describing bond-centered stripes. 
This does not imply that the stripes cannot be bond- 
centered, but rather implies that it is not sensible to de- 
scribe bond-centered stripes through spin-only models of 
coupled two-leg ladders. 

For a model of site-centered stripes described by anti- 
fcrromagnetically coupled three-leg ladders, as shown in 
Fig. |3Jd, the critical coupling required for long range or- 
der to set in is by contrast = 0. The reason is simply 
that there is no need to close a gap, as the three-leg lad- 
ders are individually gapless^^. A conventional spin wave 
analysis for such a spin-only model of three- and four-leg 
ladders was performed by Yao, Carlson, and Campbell^, 
who found that their approximation agrees reasonable 
well with the experimental data if they take J' = 0.05 J 
and J' = — 0.09J for coupled three- and four-leg ladders, 
respectively. The calculation we present in Appendix iBl 
however, singles out J' = 0.07 J for the antiferromagnctic 
coupling between spin stripes modeled by three-leg lad- 
ders. In our previous work^, we announced that a fully 
consistent spin wave theory of bond operators represent- 
ing the eight-dimensional Hilbert spaces on each rung 
of the three-leg ladders agrees perfectly with the exper- 
imental data if and only if the correct, calculated value 
J' = 0.07J is used for the coupling. 

In this context, one may ask whether it might be pos- 
sible to obtain an equally valid description in terms of 
bond-centered stripes modelled by four-leg ladders. We 
believe the answer is no, as the width of the charge stripes 
in between the ladders would be zero, and one would have 
to assume that the antiferromagnctic coupling J between 
neighboring sites of the original t-J model, would turn 
into a weak ferromagnetic coupling between the four-leg 



ladders. There would be no foundation for such an as- 
sumption. Furthermore, since the four-leg ladders are 
gapped, just as the two-leg ladders are, one would need to 
fine-tune this ferromagnetic coupling to exactly the point 
where the gap closes and magnetic order with k = 7r±7r/4 
emerges. So regardless of the agreement with the mea- 
sured spectrum one might be able to obtain, we believe 
that a spin-only model of four-leg ladders would not con- 
stitute a valid theory. 

In this paper, we provide a full and detailed account of 
our analysis of our spin-only model of coupled three-leg 
ladders. The paper is organized as follows. In section [III 
we introduce a basis as well as a set of bosonic creation 
and annihilation operators for the three-site rungs of the 
ladders, in terms of which we write both the rung Hamil- 
tonian and the spin operators on the individual sites. In 
section IIIH we couple the rungs both along the ladders 
and across neighboring ladders, and sclf-consistently de- 
termine the fiducial state such that all the terms linear 
in a single creation or annihilation operator in the re- 
sulting Hamiltonian vanish. In section IIVI we rewrite 
this Hamiltonian in terms of momentum space opera- 
tors, and expand it to bi-linear order in term of those. 
We then solve for the low energy spectrum using a multi- 
dimensional Bogoliubov transformation in section [V] In 
section IVIl we compare our results to the experimental 
data obtained by Tranquada et alr^, and investigate the 
dependence of the spectrum we obtain on the value of 
the inter-ladder coupling J'. Finally, we present our con- 
clusions in section IVlIl 



II. BASIS STATES FOR THREE-SITE RUNGS 

To begin with, consider a single rung of a three-leg 
ladder, consisting of spins which are antiferromagneti- 
cally coupled with a coupling J we set to unity (see Fig. 
H]). For later purposes, let us consider a rung belonging 

t 

Sl S2 S3 
FIG. 4: Single rung on sublattice A 

to sublattice A, i.e., set up conventions the rungs belong- 
ing to sublattice A will inherit in the following sections. 
Denoting the SU(2) spin operators for the three spin i's 
on the sites by Si,S2, and S3, the Hamiltonian for the 
rung reads 

H-^ = S1S2 + hh- (1) 

Diagonalization yields the following eigenvalues and 
eigenvectors: 

^ ^ _^ f 1^-1/2) - (im> - 2|;n> + iut» 
\\bi/2) --45(i;n>-2int> + im» 
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(2) 



Note that the two states |a_i/2) and |a_i/2) are antisym- 
metric under spacial reflections interchanging sites 1 and 
3 on the rung, while all other states are symmetric. This 
distinction will prove useful when expanding the Hamil- 
tonian for the coupled ladders in Sec. IIVI below. 

We denote the orthonormal basis formed by these eight 
states by 

M = {|^'-l/2>, |&l/2>, |a_i/2), lai/2), 

|C-^3/2>, IC-1/2), lci/2), |c3/2>}- (3) 

In this basis, the Hamiltonian matrix is trivially given by 

= -(\b-i/2){b-i/2\ + \bi/2){bi/2\ 

+ ^(|C-3/2>(C-3/2| + |C-1/2)(C_1/2| 

+ |C_1/2>(C_1/2| + |C3/2>(C3/2|). (4) 

Neither of these exact eigenstatcs, however, is suited 
as a fiducial state for spin wave theory. We are hence led 
to define a vacuum state 
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V3 
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(5) 



which interpolates between the quantum ground state 
1^-1/2) of the isolated rung with S*^ = — ^ for </) = 
and the classically Neel ordered state liti) for = 
arctan(-^) = 0.6155. The parameter </) will depend on 
the coupling between the rungs and the ladders. The 
motivation for introducing the state 1 6-1/2) will become 
clear as we determine (j) self-consistently below. 

Since we wish 1 6-1/2) to be one of our basis states, we 
replace ^ by 
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lc-3/2), 
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FIG. 5: Bosonic operator for sublattice A. 



and u = cost/), v = sin^. The relevant terms in the 
Hamiltonian ([4]) transform into 




(8) 



As a next step, we introduce bosonic creation and anni- 
hilation operators = |a_i/2)(6_i/2| etc. , as indicated 
in Fig. O The subscripts of these operators refer to the 
change in the z-component of the total spin on the rung. 
Note that these operators do not obey the commutation 
relations of independent ladder operators, as we can cre- 
ate only one "particle" with either aj or a\ or any other 
creation operator from the "vacuum" state |6_i/2). 

Completeness and orthonormality of the basis ^ im- 
plies 



|6-i/2)(6-i/2| = (l-6l6i- 
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(9) 



With and ©, the rung Hamiltonian (g]) may be 
rewritten in terms of the bosonic operators: 



1 3 . 

2-r 



-UV (^c; + Co j + - (U - V ) C^Cq 





- a\a 







lC_i + c|ci + C2C2 



(10) 



4 



5- 



5 = ^ 
antisym. 



5 = i 
sym. 



sym. 




FIG. 6: Bosonic operator for sublattice B. 



creation and annihilation operators: 







71^1 + 4,15'^ 




1 I 1-. 1 

+ «a,230l«0 ^ 






+ ^Q,84C2'^1 ' 




clc_ 


1 + '^i264co 


+ '5 0,46*^1^0 


c\cq 


+ ^0,87^2^1: 
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b\h - ala„ 





CqCq c\ci cLiC_i C2C2) 



■^ctjeeCoCo 



■^0,774^1 + 50,88^2^2 



+ Sa,13(«0 + «o) +4,16(4 + Co) 

+ S^,27(4ci + 4^1) + Sa,24(4ai + "I^l) 

+ Sa,36(aoCo + coflo) + Sa,47(4ci + c{a^) 



(14) 



The matrix elements 



On sublattice B, we introduce a similar basis Af^, with 
the only difference that the fiducial state I&1/2) has 5*^ = 
i instead of = — ^ for I6-1/2) on sublattice A: 



{|Ai)^;M = l,...,8} 

{1^1/2), |fc-i/2>,|ai/2),|a-i/2), 

IC3/2), lci/2), |c_l/2), |c_3/2>}, (11) 



with 



(\h,2) 

\\C1,2) 




(12) 



We introduce a second set of bosonic creation and anni- 
hilation operators = |ai/2)(^i/2l c^c. , as indicated in 
Fig. El 

The Hamiltonian H'^ for a single rung belonging to 
sublattice B is in analogy to (fTO]) given by 



1 3 
- + -f 

2 2 

3 



+ l^uv (Cl + Co) + - {u" - v^) ClC^ 



(13) 



For later purposes, we write the spin operators = 
± is\ and for the individual sites a = 1,2,3 on 
rungs belonging to sublattice A in terms of our bosonic 



with T = +, — , z. 



(15) 



and l/i)^ as defined in ([6]) are written out explicitly in 
Appendix [X] 

Similarly, the individual spin-operators S*^ on rungs 
belonging to sublattice B are given by: 

— 'S'q, 21-^-1 + '5'a, 41^-1 + '^0,71^-1 + '^a.li^l 
+ 82^-2-^-1 + 'S'q, 23-^-1^0 + '5'q, 43^-1^0 

+ ^a.Ti^A-^a + '^Q,84C'-2^-l + "5*0,35^0^1 

+ ^a.&h^\^\ + '5'q, 26^-1*^0 + '5'a,46^-lCo 



+ '5'a, 76^*11^^0 + '5a,87C-2C-l I 



•Sq — ('5'a ) 



(16) 
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SI 



(1 - - - Al^A_^ 



"^0^0 



+ 22-^-1-^-1 + '5'q, 55^1^1 + S^^^q^CqCq 

+ 77^-1^-1 + '^'0,88^-2^-2 

+ '^'0,13(^0 + ^0) + 16(^0 + C'q) 



■^0,27 



S. 



a,24V^-l^-l+^-lS-l) 

+ '^'0,36(^0^0 + '^d^o) 

+ '^'0,47(^-1^-1 + c-1^-1)' 
with S*^ = (mI'S'q 1'^)'' likewise given in Appendix [K[ 
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FIG. 7: The microscopic model for site centered spin stripes, 
with intra-rung couplings set to unity, inter-rung-intra-ladder 
couplings J, and inter-ladder couplings J'. The real space 
unit cell contains two rungs and is indicated by the shaded 
area in gray. 



III. COUPLING THE RUNGS 

As a microscopic model for site centered spin stripes, 
we couple the three-site rungs into three-leg ladders, 
with the spins coupled antiferromagnetically with J along 
the ladders and with J' between neighboring ladders, as 
shown in Fig. [71 The sublatticc indices assigned to each 
rung alternate in both directions, i.e., under translation 
by either of the primitive lattice vectors x = (4a, 0) or 
y = (0, a), where a is the lattice constant wc herewith 
set to unity. The microscopic model is hence given by 
the Hamiltonian 



H — ^ I Hf- + J ^ SaiSai+y + J'^ZiSi 



i-{-x 



ieA 



+ J2\^f + 'JJ2 Sajs^j+y + J' 53,5^+3, I (17) 

jeB \ a=l / 

When evaluating the spectrum of (fT7|) below, wc will set 
J = 1 (and thereby equal to the intra-rung couplings). 
The inter-ladder coupling J' is determined numerically 
by comparing the energy difference of a t-J model with a 
charge stripe with and without frustrated boundary con- 
ditions to the corresponding difference in a Hcisenberg 
model, in which the charge stripe has been replaced by 
an antiferromagnetic coupling J' (see Appendix IB|) : we 
find J' = 0.07. For the moment, however, we keep the 
inter-rung and inter-ladder couplings J and J' as free 
parameter, as this makes it easier to trace the individual 
terms in the expansion below. 

The next step is to expand p7)) in terms of our bosonic 
creation and annihilation operators, using (|14p . ()16p . and 



■ S%Sa. 



(18) 



We keep only terms up to second order in the operators. 
Since (and 5*^) contains a constant term with coeffi- 
cient II as well as the linear terms 



■Sa,i3(aS + ao) 



'a, 16 



(4 +Co), 



the expanded Hamiltonian will contain the linear term 
^ '5'q,iiSq,i6 + '^'('S'l.iiSa^ie + iq)] (cq + cq) 

Q = l 



(19) 



(20) 



in addition to the linear term 

^uv{cl+co) 

already contained in (fTU)) for each rung i on sublatticc A. 
The terms proportional to 

{al + a„) 

cancel since sf ^3 = — s| ^3 and s|^3 — 0. This cancella- 
tion can also be inferred from symmetry considerations, 
as elaborated in the following section. 

We eliminate the linear terms and (^0)1 by adjust- 
ing the parameter 0, i.e., by solving 



-uv 
2 



3 

2J E ^S,ii4.i6 + 2 J'5j,iis^,i6 = (21) 

Q = l 



with J = 1, J' = 0.07, u = coscj), v = sin^, and the 
matrix elements as given in Appendix |^ in terms 

of u and v. This yields 



V2^ 



0.5019 or 



0.3410. 



(22) 



On a formal level, the reason for rotating our basis states 
via ([5]), ([7]) and (fT2|) to begin with was that this created 
the linear term pop in H-^. Without this term, there 
would have been no way to eliminate (fT9|) . and the basis 
set would have been highly impractical for further anal- 
ysis. On a physical level, the spontaneous breakdown 
of the SU(2) spin rotation symmetry leads us to expect 
that the fiducial state of the rungs is much closer to the 
classically ordered Neel state | fit) than |&-i/2). Not 
surprisingly, the spectrum evaluated below is gapless at 
some point in the Brillouin zone, as required by Gold- 
stone's theorem for a state with a spontaneously broken 
continuous symmetry, if and only if (j) assumes the value 



IV. EXPANDING THE HAMILTONIAN 

To evaluate the spectrum of (flT)) . wc first define mo- 
mentum space operators on sublatticc A according to 
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FIG. 8: The reduced Brillouin zone corresponding to the real 
space unit cell indicated in Fig. [7] contains only l/8th of the 
full Brillouin zone and is indicated by the shaded area. 



where where N denotes the number of rungs and the 
sums over k are taken over the reduced Brillouin zone 
indicated in gray in Fig. [S] Similarly, for the creation 
and annihilation operators on sublattice B we introduce 



(24) 



S| 24 ^^'^ '^2,24 



while 



in Appendix VK\ {e.g. s\ 24 = 
Si ll = '^'l.ii) or from a symmetry consideration. As the 
Hamiltonian is invariant under reflection symmetry inter- 
changing the outer chains of each three-leg ladder (i.e., 
sites 1 and 3 on each rung) , there can only be terms con- 
taining an even number of the antisymmetric operators 
aj, flg, a|, flj^, AI]^, or A_i in the expansion. 

In addition to this reflection symmetry, we have the 
SU(2) spin rotation symmetry of the Hamiltonian (|17p . 
The spin symmetry implies that the z-component of the 
total spin, 



St 
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S. 



a J 1 



(25) 



must be conserved. This means that to second order in 
the creation and annihilation operators, only operators 
which change or S^ by the same integer can appear 
in each term. For example, we can have a term biCi or 
b\cLi, but not fo|cg. Both symmetries together imply 
that to second order, the Hamiltonian (|17|) decomposes 
into terms which contain only operators belonging to one 
particular group. 



H — Eq + HaO + HcO + Hal + H, 



c2 



Hi 



bl,cl,c-l, 



(26) 



which differ from (|23|) only in that the sign of the phases 
is reversed. Since we are only interested in the one- 
particle spin wave spectrum, we neglect the effect of the 
Hilbert space restrictions for the real space creation and 
annihilation operators (i.e., that we could create only one 
"particle" per rung) on the momentum space operators. 

As mentioned above, we only keep terms up to sec- 
ond order in the creation and annihilation operators in 
the Hamiltonian. When substituting the explicit expres- 
sions (fH|) and into pT]) . we see immediately that 
many terms yield only higher orders, while others cancel. 
To begin with, since the spin flip operators only contain 
terms of first and second order in the creation and an- 
nihilation operators and are always multiplied with an- 
other spin-flip operator, we only need to keep terms of 
first order in the expressions for s J and S~ . The expan- 
sions of the s^Sp terms are slightly more complicated, as 

and S^ contain constant terms in addition to terms 
of first and second order in the creation and annihila- 
tion operators. We have adjusted the parameter (f> such 
that the linear terms in the expansion cancel. Most of 
the quadratic terms result from multiplying the constant 
term S^^ n in the expansion of S^ with the quadratic 
terms in and multiplying n with quadratic terms 

in S^. This yields seven diagonal terms (like b\bi) and 
four off-diagonal terms (like b\ci + c\bi) for each sub- 
lattice. Three of the off-diagonal terms, those linear in 
the antisymmetric operators a\ or ai, vanish. This can 
be seen either from the explicit coefficients written out The 6x6 matrix H^ consists of the fc-independent diag- 



where Eq is a contribution to the ground state energy, 
Hao contains only the operators Oq, Oq, Aq, and Aq, and 
so on. The low energy physics we are interested in is 
contained in Hhi^ci.c-i, which we will analyze in detail 
below. As for the other terms, explicit expressions and 
expansions in terms of creation and annihilation oper- 
ators are given in Appendix [C] Hao and HcO describe 
almost dispersionless modes with energies of around 2.1 
and 2.7 (in units of Joxp which we eventually set to 
^cxp = 140 me V). Hai describes a weakly dispersing 
mode of energy of about 2.0, with a bandwidth of about 
0.2. Hc2 describes a completely dispersionless mode with 
energy 3.12. Cuts of the dispersions of these modes are 
shown in Fig.[T6]in Appendix [Cl Since these modes occur 
at energies at which we consider our spin wave theory no 
longer reliable, we will not discuss them further. 

To evaluate the spectrum of Hbi^ci,c-i, we write 



Hw,ci,c-i = i^i^k^k Jtr(iffc) 



(27) 



where 



= {B.k,bU,C_i,k,4,k,Cl-^^c_i^_^f. (28) 



7 



onal terms 



with 



'53,ll('^l,22 
3 



2 ^ 



+ 2J' S3^ii(<S'i_77 — 11), 

3 ^ 

+ 2J'sl,,iSl,,-Sl,,). (29) 
The off-diagonal terms are of the general form 

Hkaj = i^fcj-. = + cos(4fc,) + cos(fcy). 
The fc-independent coefScients 

3 

iJ?3 = - 2 J ^ <,11^^,27 + 2 J's^,n^?,27 (30) 
a=l 

result from the s^S^ terms. Expansion of the s+S*^ 
terms yields the coefficients 



TTX 


= f 


'S3,21'5'l,21' 


TTX 

^14 


= J' 


''3,71'-'l,21i 


TJX 

^16 


= J' 


*3,15'-'l,21' 


TJX 


= J' 


■^3,21 "^l^Tl' 


rjx 


= J' 


*3,21'-'l,51' 


TJX 

^3i 


= J 


*3,71''l,71' 


TJX 


= J 


*3,15'-'l,71' 


TJX 


= J 


- c+ 
*3,71''l,51' 


TJX 


= J 


*3,15'^l,15' 









and 



3 3 



a— 1 Q— 1 

^16 — J *Q,15'^Q,21i ^23~ X *q;,12'5'q,17' 



a— 1 Q— 1 

3 3 

^25^'^ X *a,21'S'a,15' -^34 = '^ X *Q,71'5'a,71 ' (^2) 
a— 1 Q— 1 

3 3 

^36^'-^ X '^a,15'5'Q,71: -^45 ~ ^ S '5a,17'5'a,51 ' 
a— 1 a— 1 

3 ^ 

All other off-diagonal elements of Hk ij vanish. 



V. SOLUTION BY BOGOLIUBOV 
TRANSFORMATION 



The Hamiltonian ([27]) can be diagonalized with a 2n di- 
mensional Bogoliubov transformation^. We begin with 
a brief review of the formalism. 

At each point in fc-space, we wish to write the Hamil- 
tonian in terms of a diagonal matrix 17, 



(33) 



^ = Mr, n = m'^hm. 



(34) 



The components of V satisfy the same commutation re- 
lations as the components of ^: 



(35) 



with 



T = diag (1,-1, 1,-1, -1,1). (36) 

This implies 



^ [fz, f J„] = J2 MuTi^M 

l,m l,m 



or 



T^MTAP. (37) 
Multiplying ^ from the right by HM yields with ((211) 
THAI = (38) 

or in components 

J2{TH)uMi,^M,,{Tn),,, (39) 

i.e., the j-th column of M is given by an eigenvector of 
TH with eigenvalue TjjD,jj. This specifies M up to the 
normalization of the eigenvectors. To obtain the normal- 
ization, it is propitious to rewrite pjp as 



(40) 



(To obtain gOl), muhiply ^ by TA/"! from the left and 
by TA/f from the right and use = 1.) Each column j 



of A'/ij must hence be normalized such that 



(41) 



Diagonalization of ()27p using this formalism at each 
point in k space with 



^ k 

yields 

Hbl.cl,c-1 



('>'l,fc' '^2,fc' '~f3,k' T4,fc' Ts^fci 7(1, fc) ' 

(7i,fc, 72,fc. 73,fc: 7lfc' 75%' 76,fc)^> (42) 

^k,i llklt,k + \ (^fc,» - ^k,i^> 



k.i 



(43) 



This Hamiltonian describes three two-fold degenerate 
modes cjk.i, which we have plotted assuming Joxp. = 
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FIG. 9: (Color online) Modes described by Hti.ci,c~i plotted 
as cuts (a) along {kx,Tr) and (b) along {n,ky) using Jcxp ~ 
140 meV. 



140 meV as cuts along (kxjTr) and (tt, fcj,) in Fig. [HI The 
two-fold degeneracy of each mode corresponds to spin 
waves with = ±1. Since we expect our spin wave the- 
ory to be reliable only for energies up to Jcxp., we will 
disregard the higher modes along with those analyzed in 
Appendix [Cl The lowest mode Wfc.i — ujk,2 ■= oj{k) is 
shown as a 3D plot for half of the reduced Brillouin zone 

in Fig. nni 

The Hamiltonian (j43p further contains a contribution 



Ebl,: 



k,ii) 



-0.22116 iV (44) 



to the ground state energy. (Here N denotes the num- 
ber of rungs. The sum extends over ^ values for k.) 
In Appendix |D1 we obtain the ground state energy by 
adding this contribution to the contributions of the other 
terms in (|26p given in Appendix [Cl and obtain Eq ~ 
— 1.73378 A^. We find that this number is in excellent 
agreement with what we would expect from diagonalizing 
the model for a finite cluster with unfrustrated boundary 
conditions. This confirms the validity of our analysis. 



VI. DISCUSSION OF THE RESULTS 

A. Agreement with the experimental data 

The significance of our results emerges in the context 
of a comparison of our spectrum with the experimen- 
tal data obtained by Tranquada et alr^ through inelas- 
tic neutron scattering on the stripe ordered compound 
Lai.875Bao.i25Cu04. The data points with their corre- 
sponding error bars are shown as black or blue crosses in 
Figs.[T^ and[T]3. respectively. 
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FIG. 10: (Color online) The dispersion ij{kx,n,y 
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In Fig. [T|d, which is directly reproduced from Tran- 
quada et ali^, the neutron data are superposed with the 
spectrum of the triplet excitation of an isotropic two-leg 
Hcisenbcrg ladder, which models bond centered stripes at 
accordingly high energies. In Fig.[T^, we have superposed 
cuts of the lowest mode u>{k) along {kx,Tr) and (tt, fcj,) 
with the experimental data. (The superposition of cuts of 
our spectrum in the x and y direction reflects the assump- 
tion that a superposition of domains with stripes along 
the two principal lattice directions has been observed in 
the experiment.) We believe it is fair to say that up to 
energies of about 180 meV, the agreement is excellent. 
(Above these energies, or more precisely above energies 
of order J, a pcrturbativc spin wave analysis becomes 
unreliable.) Likewise, the constant energy slices of the 
neutron scattering intensities x"^ ik,(-^) obtained with 
the matrix elements calculated in Appendix [E] shown in 
Fig. [11] agree very well with the experimentally mea- 
sured constant-energy slices of the magnetic scattering 
in Lai sTsBao 125CUO4 shown in Fig. 2 of Tranquada et 



B. Dependence on the inter-ladder coupling J' 

The good agreement of our results with the experimen- 
tal data up to energies even larger than J = 140 meV 
(where we would expect that the perturbative spin wave 
analysis becomes unreliable) is somewhat surprising. 
While any explanation in terms of bond-centered stripes 
through coupled two-leg laddersi^i^i^i gives immediately 
a roughly adequate estimate for the saddle-point energy 
in terms of the triplet energy gap of the individual two- 
leg ladders, it is far from obvious that a model of coupled 
three-leg ladders, which are individually gapless, should 
give a saddle-point energy consistent with the data. In 
our model, the saddle-point energy depends significantly 
on the coupling J' between the ladders. Fortunately, 
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FIG. 11: (Color online) Constant energy slices of the neu- 
tron scattering intensity (k,!^) (see Appendix [El below') 
for Jcxp ~ 140 meV and J' = 0.07 Jcxp in the magnetic Bril- 
louin zone. In the indicated energy range, only the lowest 
mode shown in Figs. 1 101 and [Tk contributes. We have replaced 
the 5-functions in frequency by Lorentzians with half-width 
A — 0.05 Jcxp and averaged over both stripe orientations {i.e., 
horizontal and vertical). 



however, it is possible to determine J' rather accurately 
through numerical comparison of a t~J model with a site 
centered spin and a site centered charge stripe to a model 
with three-leg Heisenbcrg ladders coupled by J', as de- 
scribed in Appendix |B] This analysis docs not only pro- 
vide us with the value J' sa 0.07J, but also shows that 
this value is rather robust in the sense that it does not 
significantly depend on the details of how we localize the 
stripe. To obtain a better understanding of the depen- 
dence of our final results on this coupling, we have ob- 
tained the spectrum for a number of different values of 
J' by solving (|2ip numerically for each value, and pro- 
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FIG. 12: The saddle point energy lj{-k,-k) of the lowest mag- 
netic mode uj{k) calculated for various values of J' / J (points) . 
Fitting yields u}{tt,tt) 1.47V J' J (solid line). 
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FIG. 13: Auxiliary models of (a) two weakly coupled three- 
leg ladders and (b) two weakly coupled spin i chains used in 
the discussion to understand the square root dependence of 
aj(7r,7r) on J' depicted in Fig. 1121 



cccding with the Bogoliubov transformation with the re- 
sulting values for u(J') and v{J'). The results for the 
saddle-point energies w(7r, tt) are shown in Fig. [T^ (black 
dots). Fitting the data yields 



uj^n, tt) 



(45) 



to an excellent approximation up to values where J' be- 
comes comparable to J (solid line). 

The square root dependence of a;(7r, tt) on J' can be un- 
derstood by considering a model of two three-leg ladders, 
which are weakly coupled by J', as shown in Fig. 113b ,. 
The low energy excitations of the individual three-leg lad- 
ders are spin i spinous, which are gapless. The coupling 
J' induces a linear confinement potential 



Viy) - F\y\ 



(46) 



between pairs of spinous, since the links coupling the 
chains effectively become decorrelated in the region be- 
tween them. The situation here is similar to a system 



10 



of two coupled spin ^ chains shown in Fig. [T3b . where a 
weak coupling Jj^ between the chains is known to induce 
a linear confinement potential between pairs of spinons^^. 
In the model of coupled chains, the force between the 
spinous is proportional toi^i^ 

F w {SS)±Ji_ cx Jl/J. (47) 

For the model of two coupled thrcc-lcg ladders we con- 
sider here in the context of understanding the dependence 
(|45p of our spin wave analysis, however, we assume that 
the spin correlation (SS) between the sites coupled by 
J' is to lowest order independent of J'/ J. Therewith 
we account for the static correlations present due to the 
long range order we assume. For the confinement force 
in our auxiliary model of two coupled three-leg ladders 
we hence assume 

F cx J'. (48) 

The spinon confinement will then induce a gap A, which 
corresponds to the the ground state or zero-point energy 
of the linear potential oscillator for the relative motion 
of the spinous. The dispersion of the spinous is linear 
for both the individual spin chains and the individual 
three-leg ladders^, 

eiky) « v\kyl (49) 

with w oc J as J is the only energy scale is these models. 
(In (|49p . we have shifted both spinon branches to the 
origin.) The ground state energy Eq of a constant force F 
oscillator of linearly dispersing particles with velocity v, 
however, is proportional to V Fv—. This implies A cx Jj_ 
for the two weakly coupled chains and A cx ^/J'J for the 
two weakly coupled three-leg ladders with the additional 
assumption of static correlations sue to long range order. 

To see why this gap A corresponds to the saddle point 
energy w(7r,7r) in the spin wave analysis above, con- 
sider the transformation properties of our auxiliary model 
of the two weakly coupled three-leg ladders shown in 
Fig. 113b . under the parity reflection x — > —x. The gapped 
spinon-spinon bound state is odd under this symmetry, 
which in the language of momenta of the site centered 
stripe model corresponds to a shift of j . Since the ground 
state of the stripe model has order with k^ = tt ± ^ , the 
gapped excitation will correspond to k^ = tt. Following 
this line of reasoning, we can understand the square root 
dependence of ^(Tr, tt) depicted in Fig. [T^l 

VII. CONCLUSIONS 

In this work, we have provided a full and detailed ac- 
count of a spin wave analysis of a coupled three-leg ladder 
model for spin stripes in copper oxide superconductors. 
We have numerically evaluated the inter-ladder coupling 
J' induced by the charge stripes in between both site- 
centered spin stripes modelled by three-leg ladders and 



bond-centered spin stripes modelled by two-leg ladders. 
As reported previously^, for the latter we obtain a fer- 
romagnetic coupling J' = —0.05 J, which is not sufficient 
to close the energy gap of the individual two-leg ladders. 
This does not imply that the stripes cannot be bond- 
centered. It docs imply, however, that a description in 
terms of spin-only models of coupled two-leg ladders is 
not sensible. 

For site-centered spin stripes modelled by three- 
leg ladders, we obtain an antiferromagnetic coupling 
J' = 0.07J. We have calculated the spectrum, the 
staggered magnetization, the dynamical structure factor 
(fc,a;), and the ground state energy for a spin-only 
model of coupled three-leg ladders using a linear spin 
wave analysis of bosonic operators representing the 
eight-dimensional Hilbert spaces on each three-site 
rung. The analysis makes no assumptions except for the 
model itself, and contains no variational parameter, as 
even the inter-ladder coupling J' is evaluated through 
exact diagonalizations of small t-J clusters with and 
without charge stripes in between the spin stripes. We 
find excellent agreement with the experimental data of 
Tranquada et al^. The experimental data hence point 
towards site-centered, and not, as previously asserted, 
bond-centeredi^i^^i^ stripes. 
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Appendix A: Matrix elements of tlie individual spin 
operators on rungs 

In Sec. ini we have written out the spin operators 
and with t = +, — , z on the individual sites a = 1,2,3 
on rungs belonging to sublattice A and B, respectively, 

in the basis sets and specified in ([6]) and (llip . 
The matrix elements 

sl,^. = (A3) 
SI,. = {l^\SlWf. (A4) 

are explicitly given by: 
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Appendix B: Estimation of the coupling J' between 
ladders 



The effective, antiferroniagnetic coupling J' between 
neighboring spin stripes described by three-leg ladders is 
induced by the charge stripes between the ladders, since 
an antiferromagnetic alignment of the spins on each side 
of the charge stripe enhances the mobility of the holes. 
To determine this coupling, we have exactly diagonalized 
16 site clusters of itinerant spin 1/2 antiferromagnets de- 
scribed by the t-J model^i^ with J = OAt, two holes, and 
periodic boundary conditions (PBCs), in which the spin 
stripes are localized. The Hamiltonian is given by 



t-J-B = -t Y.VG{cl^C^,,+\,.C.)Fo + JY.S^S, 

zGshadcd area 

where the first two sums extend over all nearest-neighbor 
pairs and Bj = ±B denotes a staggered magnetic 

field, as indicated by the signs in Fig. [141 The Gutzwiller 
projector Pg eliminates doubly occupied sites. 

We compare the ground state energies we obtain for 
clusters with unfrustrated PBCs shown in Fig. [Tik with 
the ground state energies we obtain for clusters with frus- 
trated PBCs shown in Fig. [T4b . in which the 16-site unit 
cells on the right are shifted by one lattice spacing to the 
top. We then consider spin-only Heisenberg models of 
three-leg ladders, which consist of only the sites in the 
shaded areas in Figs. [T4k and [T4b . subject to the same 
staggered field Bi, and couple them antiferromagnetically 
by J', as indicated. The Heisenberg models are described 
by 
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FIG. 14: Finite size geometries with (a) unfrustrated and 
(b) frustrated periodic boundary conditions for site-centered 
stripe models. The spin stripes are localized by a staggered 
magnetic field B, as indicated by the signs in the grey shaded 
areas. 
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where Jij = J for all nearest-neighbor links inside the 
shaded areas, but Jij = J' for nearest-neighbor links 
across the horizontal boundary lines between those areas. 
We again compare the ground state energies for unfrus- 
trated PBCs, where J' couples sites 10 and 1, 11 and 2, 
etc. for the three-leg ladders shown in the shaded areas 
in Fig. [Hk . with frustrated PBCs, where J' couples sites 



11 and 1, 12 and 2, etc. , as shown in Fig. [Mb. Finally, 
we determine J' such that the difference in the ground 
state energies between frustrated and unfrustrated PBCs 
in the t-J clusters matches this difference in the spin-only 
ladder models. 

With B = 0.170J we obtain J' = 0.071J, as detailed 
in Tab.lll The value for B is chosen self-consistcntly such 
that the mean value of the staggered magnetization we 
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Fig. [Hi 


Fig.IHb 




t~J model with spin and charge stripes 


t = 2.5, J = 1, B = 0.17, iV = 16, 2 holes 


Et-J-B 


-21.3409 


-21.2405 
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Eb 


-0.5415 






{SI) 


-0.3140 






(SI) 


0.2411 
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0.9437 






(ns) 


0.8781 






three-le 


g Heisenberg ladder with coupling J' 


J' = 


= 0.071, J = 1, B = 0.17, iV = 


12 


Ej-B 


-8.0068 


-7.9065 


0.1003 



TABLE I: Ground state energies and expectation values of 
the staggered magnetic field term Hb = "^BiSf, of St, and 
of the electron densities fn obtained by exactly diagonaliza- 
tion of the t-J model (IBlfl as well as the coupled Heisenberg 
ladders described by (|B2|) for the clusters shown in Fig. 114b , 
and 114b for unfrustrated and frustrated periodic boundary 
conditions, respectively. 



obtain for our spin wave theory in Appendix |F1 



2 (-(^0 + (^2)) = 0-2903, 



(B3) 



matches the corresponding value in the t-J cluster shown 
in Fig. 114b .. For this cluster, however, there are two values 
for the staggered magnetization, depending on whether 
we consider the overall magnetization 



\{-{si) + {si) 



0.2775, 



(B4) 



or the magnetization on only those sites which are occu- 
pied by electrons (which differs since the holes are not 
strictly localized on the chain in between the shaded ar- 
eas in Fig. \Wd]: 



("1) ("5), 



0.3037. 



(B5) 



We assert that these two values constitute lower and up- 
per bounds of what we would expect in a spin-only model, 
and chose B such that the spin wave theory gives the 
mean value of these bounds (this value is 0.2906). 

For completeness, we also provide the details of the 
corresponding calculation for the effective coupling J' 
between bond-centered stripes modelled by two-leg lad- 
ders, in Tab [III The calculation differs only in that we 
now use the finite size clusters depicted in Figs. [T5k 
and [TSb . The value of the staggered magnetic field 
B ~ 0.225 is chosen such that the staggered magnetic 
field energy Eb described by the last term in (|Bip is equal 
to the value we obtained for the three-leg ladder model 
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FIG. 15; Finite size geometries with (a) unfrustrated and 
(b) frustrated periodic boundary conditions for bond-centered 
stripe models. The spin stripes are localized by a staggered 
magnetic field B as indicated by the signs in the grey shaded 
areas. 



listed in see Tab. HI We obtain a ferromagnetic coupling 
J' = — 0.051J. This value is not large enough to close 
the spin gap A « J/2 oi the the individual two-leg lad- 
ders, and hence precludes a description of bond-centered 
stripes in terms of spin-only models of coupled two-leg 
ladders^^. 

It should be born in mind that the values we obtain 
here are only estimates, as it is impossible to calculate a 
precise value for a coupling between spin-only models of 
stripes, as these constitute a rather crude approximation 
themselves. We are confident, however, that the coupling 
for three-leg ladders modelling site-centered stripes is be- 
tween J' = 0.05J and J' = O.IJ, and that the absolute 
value of the coupling between two-leg ladders describing 
bond-centered stripes is significantly smaller than this 
coupling. 



Appendix C: Analysis of Eo, Hao, Hco, Hai, and Hc-i 

In this appendix, we expand and analyze the terms in 
the Hamiltonian (|26p which have no influence on the low 
energy spectrum in this appendix. These terms, however. 
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Bond-centered stripe models 


PBCs 


unfrustrated 


frustrated 


A 


as in 


Fig. [Hi 


Fig.IHb 




t~J model with spin and charge stripes 
t = 2.5, J = l,B = 0.225, N = 1&,2 holes 


Et-J-B 

Eb 


-21.3428 
-0.5395 


-21.2526 


0.0902 


two-le^ 
,7' = 


; Heisenberg ladder with coupling J' 
-0.051, J = 1, B = 0.225, = 8 


Ej-B 


-5.1557 


-5.0655 


0.0902 



TABLE II: Ground state energies and the expectation value 
of the staggered magnetic field term Hb = X^-BiSf obtained 
by exactly diagonalization of the t-J model (|B1|I as well as the 
coupled Heisenberg ladders described by (|B2|I for the clusters 
shown in Fig. 115b and 115b for unfrustrated and frustrated 
periodic boundary conditions, respectively. 



are required for evaluations of the ground state energy. 

1. The bare ground state energy Eo 

The bare ground state energy i?o accounts for the con- 
stant terms in the Hamiltonian p7|) . It is given by 
^0 = Y.k 2^0 = Nia with 



y 3 ^ 



a=l 



-1.45842. 



(CI) 



2. Evaluation of the spectrum of Hao 

The term Hao is given by 



+ U.k{alk + a^,k){Al^ + A^^^) , (C2) 



with 



1 3 ^ 

(C3) 

3 

^aO.fc = J'cOs(4fc:r)s3^i3S'f43-h JcOs/Cj^^S^ igS'^ ig. 

a=l 

(C4) 

To a reasonable approximation, we obtain low energy 
modes described by Hao by diagonalizing (jC2[) at each 
point in k space in the reduced Hilbert space spanned by 
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FIG. 16: (Color online) Modes described by HaO (blue), Hco 
(green), Hai (red), Hc-2 (black) plotted as cuts (a) along 
{kx,n) and (b) along (-Kjky) using Jcxp ~ 140 meV. 



where 



io) ^ n 1^-1/2). ■ n 1^1/2) 



(C5) 



is the bare vacuum unrenormalized by spin wave theory. 
This yields two almost dispersionless modes 

(Wa0,fc)i/2 = \/ ^aO + CaO.fc ± C<iO,fc, (C6) 

with energies of about Sao = 2.07, or 290 meV if we as- 
sume Jcxp. = 140 meV. Cuts of the dispersions of these 
two modes are shown in blue (color online) in Fig. 1161 
HaQ also gives rise to a contribution 



-0.00008 iV (C7) 



to the ground state energy. (Here N denotes the number 
of rungs, which implies that the sum extends over ^ 
values for k.) 



3. Evaluation of the spectrum of Hco 

A similar analysis of 



|0), 4|0), 4|6), and alAl\0) 



t 4tin\ 
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with 



a=l 

+ 2J'4ii(5f66-5fii), (C9) 



with energies of about Sco = 2.71, or 380 mcV, which is 
shown in green (color online) in Fig. [161 Hco also gives 
rise to a contribution 



Eco = J2 (^-0 - V^'o + e'cfc ) - -0.00048 N (C12) 



to the ground state energy. 

4. Evaluation of the spectrum of Hai 

The term 

Hal = E 

+ UAaik^Uk + . (C13) 

with Eai = ^ao as given in (jC3[) and 

3 

^al.fc = J'c0s(4fc2:)s;^4i5'];^41 + JcOSfcy S^Al^aAl' 

a=l 

(C14) 

can be diagonalized by a Bogoliubov transformation. We 
obtain 

(C15) 



with 



It yields a two-fold degenerate, weakly dispersing mode 
with an energy of about 1.95, or 273 meV, which is shown 
in red (color online) in Fig. I16i as well as a contribution 

Eai = J2 {^-^■'^ - ^-i) = -0-05363 N (C17) 



5. Evaluation of the spectrum of Hc2 

Finally, 

Hc2 = E^^2(4,fcC2,^ + Cl2,feC_2,J (C18) 



^co,k = J'cos{4:k^)sliQSliQ + JcoskyJ2sa,i6Sa.W' with 

(CIO) 

yields two additional, almost dispersionless modes 

(Wc0,/=)i/2 = \/ + icO.k ± CcO.fc, (Cll) 



Wc2 = ^w' + 2jE<n('5o,88-^a,ll) 



2"^ *3,ll('^l,88 "^1,11)1 



(C19) 



describes a two-fold degenerate, completely dispersion- 
less mode with an energy of uic2 = 3.12, or 436 meV. 



Appendix D: Ground state energy 

To evaluate the ground state energy Eq, we collect the 
contributions from ([CT]), jCTl, (p^ . (pJll . and (??). 
This yields 



-Bo = -1.73378 7V, 



(Dl) 



where N is the number of rungs. This number is in good 
agreement with what we would expect from the results 
of exact diagonalizations of 



(D2) 



for small clusters of 12, 18 and 24 sites with unfrustrated 
boundary conditions, as shown for = 4 rungs by the 
shaded areas in Fig. [T4k . Specifically, we obtain Eq = 
-1.827iV for a cluster with = 8 rungs, J' = 0.07, and 
unfrustrated boundary conditions with a shift of 3 lattice 
spacings. If we compare the nearest-neighbor spin-spin 
correlation we obtain from exactly diagonalizing the same 
24 site cluster with J' = 1, (S'iS'j) = -0.343, to the the 
value predicted by standard two-dimensional linear spin 
wave theory^, (SiSj) = -0.329, we are led to estimate 
that the bond operator spin wave theory developed here 
should yield a number around 

—0 329 

£:o = -1-827 iV '■ = -1.751 TV. (D3) 

-0.3432 ^ ^ 



(C16) This differs only by 1% from ([DT 



to the ground state energy. 



Appendix E: Matrix Elements 

The dynamical structure factor measured in neutron 
scattering is given by 

X+"(fe,a.) = E I (0| S+ \n) I'Siu - c.„), (El) 
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where |0) is the ground state and the sum extends over 
all excited states |n) with energy and 



-ikri ^+ 



(E2) 



is the Fourier transform of the spin raising operator Si 
at lattice site I with respect to original lattice, i.e., the 
sum runs over all lattice sites. This implies 



(E3) 




(E4) 



with Ri and Rj as indicated in Fig. [7] for sublattices A 
and B, respectively. In analogy to Fourier transforms of 
the bosonic creation and annihilation operators p3p and 
([24|) , we further Fourier transforms of the spin operators 
with respect to the rung sublattices A and B according 
to 



N ^ ' 



-ikRi -+ 



St 



(E5) 
(E6) 



and express the operator (jE2p in terms of them: 



Si 



EE^' 



^■"Sl 



E 

a 

E 



(E7) 



As we are interested only in the contribution of the low- 
energy mode Wfe^i to (^j the only matrix element 
we need to evaluate is 



14+171,01= |(0|5+7U|0>| 



(E8) 



for all values of fc. (The second low-energy mode 01^,2, 
which is degenerate with the first, does not contribute 
to (|E8p and hence to x"' (fc,i^), but instead yields a 
contribution to x ^ [k, w) which is identical to the one we 
calculate below.) Keeping only terms which contribute 
to this mode and are linear in the expansion (|14p of s+, 
we obtain for sublattice A 



or in Fourier space with p]5| . (fSSj) . and (128 



'^Q,2l'^l,fc 



„+ J 



+ S 



Q, 15^-1. -fc 



(E9) 








FIG. 17; (Color online) Numerical evaluation of the ma- 
trix elements K0|S'^|7i,fe)|^ for the entire Brillouin zone k £ 
[0, 27r] X [0, 27r]. Note the strong enhancement aground the 
antiferromagnetic ordering wave vectors k = (tt ± 7r/4, vr). 



and similarly with ((16)) and ((241) ^ 

Sk,a = SaS2^-l,k + Sa.n'^-IM + -^isi^l -fc 

= {Sli2,0,Sl„,0,Sl,i,0)-^k. (ElO) 

We then use ^Pfc = ./l/fcFfc and ((42l) to express s'^ ^ and 

•^fc.a terms of 7i,fe and jj^^, recall S'+^j = s+ j^-, and 
finally obtain 

|(0l4l7i.fc)|' = I E^"^'"''""'H<2i(^4,ii+A4,2i) 

a = l 

+ 4,71 (^^fc,31 + Mk,41 ) + S+ 15 ( A/fc,51 + Mfc,61 ) } 

(Ell) 

Note that x^^(fc,w) = as there is no term linear in b\, 
c\, or c_i in the expansion ((T4l) for s^. 



Appendix F: Staggered magnetization 

The staggered magnetizations on the outer and inner 
chains of our three-leg ladders, (sf ) = (s|) = — (S'l) = 
— (S"!) and ^Sj) = ^{^2)1 respectively, are given in by 
the the bare values n = — n with a = 1,2,3 minus 
corrections from the individual terms in the Hamiltonian 
pel) . From ((Ml) - ((XTOl) with ([22]), we obtain for the bare 
values 



('-'2) bare 



'1,11 



'2,11 



= -0.4633 
= 0.4265 



(Fl) 
(F2) 



We expect that the largest corrections arise from 
Hbi.ci,c-i, as this part contains the only low energy mode 
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of the theory. The for this part relevant terms in the ex- 
pansion (fT4|) for are 

b\bi (4.22 - 4,11) + 4ci (4,77 - 4,11) 

+ (4,55 - 4,11) + (^iCl + c\b^) 4,27 • 

For convenience, we define 

Pa,l = (■Sq,22 ~ 



Pa,2 = (4,77-4,11) (F3) 
Pa,3 — (Sa,55 ~ ^q,1i) • 

and express the operators bi,ci,c^_i as weh as their her- 
mitian conjugates via 5* = MT through the 7 operators 
(see (ESI)), (IMl), and dUl) above). Using 7^-10) = (0|7j = 
and (0|7i7j|0) = Sij, we obtain for the corrections from 

i^bl,cl,c-l 



"'-tfbl,cl.c-l ]\[ 



^ I [pa,lMfc,2lMM2 + Pa,2Mfc,4lMM4 + Pa,3Mfc,62M^ 26 + 4,27 (Mfc,2lMM4 + MkAlMl,^) 
|^PQ,lA/fc, 23^^32 + PQ,2M;e,43-/V^fc^34 + PQ,3Mfc,64A^fe^46 + Sa^27 (-^fc,23A/fc^34 + ^-^^,43^/^32) 
Pa,lMfc,26Mj62 + Pa,2Mfc,46M^64 + Pa.sMkSbMl^g + 27 (^^,26^^64 + Mk,Af>Mlf^^) 



In this sum, the terms in the first pair of square brackets 
originate from the low energy eigenmodes 71 , 72 in the 
energy range from to about 190 meV in Fig. [SI the 
terms in the second pair originate form 73 , 74 at about 
290 meV, and the terms in the third pair originate from 
75,76 at about 460 meV. Evaluation yields 

isDf, = 0.1752 + 0.0018 + 0.0003 

= 0.1773 (F4) 
(sDf, = -0.0864- 0.0210 - 0.0001 

\ -Hbl,cl,e-1 

= -0.1075. (F5) 



where 



ao,k 



[SaM *Q,llj 



CaO.fc 



-^2 



(FIO) 



(Fll) 



(F12) 



for the outer and inner chains, respectively. As expected, 
the low energy mode we compare to the experimeniji^ 
yields the dominant contribution. 

The corrections arising from Hai arc evaluated in 
complete analogy. The Bogoliubov transformation 
^'a = MaTa With 'i' a = iA-i,0-\)^ is uow Only two- 
dimensional, and the only contribution comes from the 
term 



at a, in (fT4 



^Q,ll"l"'l 



with eao and S,a0.k as given in (jC3|l and ()C4p and eco and 
^co.k as given in ((C9l) and ([CTU)) . Evaluation yields 



= 0.0001, (F13) 
(4)i?„„ + (s2^)ij,o = -0.0001, (F14) 



■^"^-Hal ^ ■Sa,ll-^-^a,fc,21^^Xfc, 



12- 



Evaluation yields 



(F6) negligible small contributions. 



Summing up (|FTI), dEl, (|F7|, (IfBI) . and (|F2|), (|F5|1 . 
F8|, (|F14p . we obtain 



= 0.0127, 
= -0.0117. 



Finally, Hao and 77co give rise to corrections 

(-4,11) TjJ^ljr 



(F7) 
(F8) 



(F9) 



(4) 



-0.2732, 
0.3072 



(F15) 
(F16) 



for the staggered magnetizations on the outer and inner 
chains of the three-leg ladders, respectively. 



oO,fc 
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